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ABSTRACT 

Context. The BEER algorithm searches stellar light curves for the BEaming, Ellipsoidal, and Reflection photometric modulations that 
are caused by a short-period companion. These three effects are typically of very low amplitude and can mainly be detected in light 
curves from space-based photometers. Unlike eclipsing binaries, these effects are not limited to edge-on inclinations. 

Aims. Applying the algorithm to wide-field photometric surveys such as CoRoT and Kepler offers an opportunity to better understand 
the statistical properties of short-period binaries. It also widens the window for detecting intrinsically rare systems, such as short- 
period brown-dwarf and massive-planetary companions to main-sequence stars. 

Methods. Applying the search to the first five long-run center CoRoT fields, we identified 481 non-eclipsing candidates with pe¬ 
riodic flux amplitudes of 0.5-87 mmag. Optimizing the Anglo-Australian-Telescope pointing coordinates and the AAOmega fiber- 
allocations with dedicated softwares, we acquired six spectra for 231 candidates and seven spectra for another 50 candidates in a 
seven-night campaign. Analysis of the red-arm AAOmega spectra, which covered the range of 8342-8842 A, yielded a radial-velocity 
precision of ~ 1 km/s. Spectra containing lines of more than one star were analyzed with the two-dimensional correlation algorithm 
TODCOR. 

Results. The measured radial velocities confirmed the binarity of seventy of the BEER candidates—45 single-line binaries, 18 double¬ 
line binaries, and 7 diluted binaries. We show that red giants introduce a major source of false candidates and demonstrate a way to 
improve BEER’s performance in extracting higher fidelity samples from future searches of CoRoT light curves. The periods of the 
confirmed binaries span a range of 0.3-10 days and show a rise in the number of binaries per Alogf toward longer periods. The esti¬ 
mated mass ratios of the double-line binaries and the mass ratios assigned to the single-line binaries, assuming an isotropic inclination 
distribution, span a range of 0.03-1. On the low-mass end, we have detected two brown-dwarf candidates on a ~ 1 day period orbit. 
Conclusions. This is the first time non-eclipsing beaming binaries are detected in CoRoT data, and we estimate that ~ 300 such 
binaries can be detected in the CoRoT long-run light curves. 

Key words, binaries: spectroscopic - binaries: eclipsing - brown dwarfs - Techniques: photometric - Techniques: spectroscopic - 
Techniques: radial velocities 


1. Introduction 


The space-based p hotometric surveys CoRoT (Rouan et al J 
1 1998 1 IBaglinl 1200 3) and Kepler dBorucki et al.ll2Qloli ~ were de¬ 
signed mainly to detect minute exoplanetary transits. Each of 
these missions has provided over 160,000 continuous stel¬ 
lar light curves with time-spans of tens to hundreds of days 
and a photometric precision of 10 3 —1 0 4 per measurement 
(lAuvergne et al.l 120091: iKoch et al.ll2010l) . Hundreds of transit¬ 
ing planets were indeed detected (e.g., iMoutou et al.l 1 201 3l 
[Rowe et all [20141) . but the unprecedented quality of the light 
curves also enabled the detection of other flux variations of as- 
trophysical origins. One such variation is the relativistic beaming 
effect. 


iRvbicki & Lightmanl ( 1979 1 showed that several factors con¬ 
tribute to the beaming effect, which increase (decrease) the ob¬ 
served brightness of any approaching (receding) light source by 
~ |4v/e/c|, where vr is the radial velocity (RV) of the source and 
c is the speed of light. Thus, for short-period (1-10 days) brown- 
dwarfs (BDs) or massive planetary companions (M p sin i > 
5 Mj) of solar-like stars, the beaming amplitudes are in the range 


of 10~ 4 -10~ 5 . lLoeb & Gaudil d2003l) have predicted that Kepler's 
photometric precision would be sufficient to detect such com¬ 
panions. For stellar binaries with typical RV semi-amplitudes 
of 10 - 100kms _1 , the be aming ampl itudes are in the range of 
10 3 —10 For this reason. IZucker, Mazeh, & Alexanderl (120071) 
have predicted that CoRoT and Kepler will also detect hundreds 
of non-eclipsing binaries of this type, and create a new observa¬ 
tional category: beaming binaries. 


Soon after the first Kepler light curves became available, sev¬ 
eral studies measured the beaming e ffect of a few echpsing bi¬ 
naries (EBs) detec ted by Kepler (e.g. . Ivan Kerkwiik et al. 20101 
ICarter et~aD 1201 lh. However, as mentioned by Loeb & Gaudil 
(l2003h and by Zucker. Mazeh. & Alexanderl (l2007l) . for short- 
period binaries and planets, the amplitude of the beaming ef¬ 
fect might be comparable to (or even smaller than) another two 
photometric effects—the ellipsoidal and the reflection light vari¬ 
ations. The ellipsoidal variation is caused by tidal interactions 
betwee n the two components of the binary (e.g., iMorri sill 98.4 
iMazehl 120081 ). The reflection variation is caused by the bright¬ 
ness difference betwee n the ’day’ side and the ’night’ sid e of 
each component (e.g.. IWilsonl 1 1 990l: lHarrison et alJl2003h . By 
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accounting for the three effects, several studies succeeded to de¬ 
tect the weak beaming effect caused by a transiting BD or even 
a tra nsiting massive planet in CoRoT and Kepler light curves 


(e.g., Mazeh & Faigler 2010k Shnorer et al. 

201l|; |Mazeh et al.l 

201211 Jackson et all20ll’ Mislis et al. 2012 

). 


To find non-eclip sing short-period beaming binaries, 
iFaigler & Mazeh: ( 201 1) introduced the BEER algorithm, which 
searches light curves for a combination of the three photomet¬ 
ric effects caused by a short-period companion—the BEaming, 
Ellipsoidal, and Reflection periodic modulations. BEER approx¬ 
imates each of the three effects by a sine/cosine function relative 
to phase zero taken at the time of conjunction—when the lighter 
component is in front of the heavier one. The reflection and the 
beaming effects can then be approximated by cosine and sine 
functions with the orbital period, respectively, while the ellip¬ 
soidal effect can be approximated by a cosine function with half 
the orbital period. 

Detection of BEER-like modulations in a stellar light curve 
does not yet prove the binary nature of the star, since sinu¬ 
soidal flux modulations cou ld be produced by other effects as 
well (e.g jAigrain et al.l2004l) . To confirm BE ER detections, one 
needs to perform RV follow-up observations dFaigler & Mazehl 
1201 ll) . The first RV co nfirmation of non-e clipsing beaming bi¬ 
naries was reported by IFaigler et al.l (120121) . Candidate binaries 
for that study were detected with BEER in Q0-Q2 Kepler light 
curves, and seven of them were confirmed using RV measure¬ 
ments. __ 

Paper I of the current series dFaigler et all2013h reported the 
discovery of Kepler -76b—the first ho t Jupit er discovered with 
BEER. In Paper II, IFaigler & Mazehl (1201 5l) showed evidence 
for equatorial superrotation of three hot Jupiters measured by 
Kepler —KOI-13, HAT-P-7, and Kepler-76b. 

In this paper we present RV confirmation of seventy new 
beaming binaries found by BEER in CoRoT light curves. The 
targets were selected from the first five long-run center CoRoT 
fields and were confirmed usin g the AAOmega multi-object 
spectrograph dLewis et al.ll2002i) . Section [2] presents the BEER 
search applied to the CoRoT light curves. Sect. [2 describes the 
spectroscopic follow-up observations, Sect.[4]details the spectral 
analysis and derivation of RVs from the spectra. Sect. [^explains 
the orbital solutions performed and the statistical methods ap¬ 
plied to separate true BEER binaries from false detections, Sect. 
[6] discusses the use of the findings to evaluate the BEER algo¬ 
rithm performance, Sect.Qfocuses on the mass ratio and orbital- 
period distributions of the new binaries, and Sect. [8] summarizes 
the findings. 

2. BEER photometric search for binaries 

To detect beaming-binary candidates, we analyzed the ~ 40,000 
white light curves of the CoRoT fields LRcOl, LRc02, LRc03, 
LRc04, and LRc05. We did not use the red, green, or blue 
light curves of targets having chro matic information (e.g., 

I Ai grain et al.l2008tlDeleuil et al .120 j~Ti) due to their lower signal- 
to-noise ratio (S/N). 

The light curve analy sis co nsisted of several steps. First, 
oversampled light curves (ISurace et al.ll2008l) were rebinned to 
512 s. Then, we corrected for jumps in all light curves. Jumps 
were identified by calculating a median filter to each light curve 
and detecting outliers in its derivative. The correction was made 
by subtracting the median filter from the light curve around 
the identified jump epoch. A cosine-transform-based detrending 
and 4cr outliers removal w ere then performed using ROBUSTFIT 
dHolland & Welsciilll977h . Next, a fast Fourier-transform (FFT) 


-based power spectrum (PS) was calculated for each light curve, 
and the fiv e most p rominent peaks were identified and analyzed 
following lFaigler et al.l (120131) . As detailed there, for each of the 
five peaks we derived the BEER amplitudes and the estimated 
mass and albedo of the presumed companion by fitting the am¬ 
plitudes with a BEER model, assuming that the orbit is circular 
and the peak corresponds to either the orbital period or half the 
orbital period. Hence, for each light curve we evaluated ten pos¬ 
sible orbital periods. 

After fitting a circular BEER model to each light curve at 
each of its suspected periods, we assigned each fit a score in 
the 0-1 range, with 1 being the best. The score of each fit was 
calculated as 


S Total - S BEER_S/nS min_ S/N S' x iS sin iS Albedo > 


( 1 ) 


where S bkkr_s/n is the BEER model S/N score, S mm _s/N is the 
score of the minimum S/N of the two BEER harmonics, S y i is 
the fit score, .S' sm , is the score of the model-derived sin;, and 
S albedo is the score of the model-derived geometric albedo. Each 
of these partial scores is the result of a dedicated scoring function 
that gives a score in the range of 0-1 for its associated parameter. 
The scoring functions we used were 


Sbeer_s/n - 1 - exp 


BEERuS/N 

Cbeer_s/n 


( 2 ) 


S min_S/N — 1 exp 


minS /N 

f inln_S/N 



(3) 

(4) 


S sin i — 1 exp 


sin; 
r 

'-'.sin i 


and 


(5) 


/ max(Albedo - C a i b _ m , 0)\ 

s Albedo = exp--- , (6) 

\ ^alb_s / 

where C B eer _S/N 5 Cmin_ s/N) C^. 2 ’ fsin i, T'aib.m, and Caib_s lire con¬ 
stants that calibrate the behavior of the scoring functions. The 
most likely orbital period of each light curve was selected as the 
period with the highest total score S Total out of its ten evaluated 
periods. 

As an illustration. Figs. [I] and [2] show the light curve analy¬ 
sis of CoRoT 105962436 and CoRoT 104674562, respectively. 
We present 70-day-long parts of the original and detrended light 
curves, the entire light curve power spectra, and the entire phase- 
folded and binned light curves, with the best-fit BEER models 
superposed. We note that the most prominent peak in Fig. Q] at 
1.11 day~' is caused by the ellipsoidal effect, whose frequency 
is twice the orbital frequency (i.e., the predicted orbital period is 
1.8 days). Figure [2] shows that the BEER modulation of CoRoT 
104674562 at a period of 4.6 days is almost buried in the noise, 
and from the PS it seems that the modulation at ~ 20.4 days is 
more evident. Nevertheless, thanks to its period-selection pro¬ 
cess, BEER found the correct orbital period in this case as well 
(as we show in Sect. [3. 

After the score assignment and the best-period selection, we 
visually inspected the 200 highest score candidates of each field 
and assigned priorities of 1-3 to targets that we deemed viable 
binary candidates, with priority 1 assigned to the best candidates. 
These priorities were assigned through visual inspection of the 
photometric power spectrum, the goodness-of-fit of the BEER 
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HJD - 51545 



x 10 - 3 105962436 



Fig. 1. BEER light curve analysis of CoRoT 105962436 (a) A 70-day long part of the CoRoT white light curve, normalized by 
its median. For clarity, the original light curve was shifted upward by 0.03 relative to the cleaned and detrended one. The cleaned 
light curve demonstrates the jumps correction, outlier removal, and detrending functionality, (b) FFT-based power spectrum of the 
detrended light curve, with its maximum value normalized to 10. Vertical dashed lines mark the first three harmonics of the candidate 
orbital frequency, (c) Phase-folded and binned light curve (blue) and the best-fit BEER models assuming a circular orbit (green). 
The residuals were shifted downward for clarity. 


model, the correlation structure of the time-domain residuals, 
and other target specific features. Naturally, the visual inspec¬ 
tion method and the resulting target priorities are subjective and 
prone to human bias and errors, but we nevertheless used it for 
lack of a better software-based method at this stage. 

In this way, we selected a total of 481 candidates for RV 
follow-up from all five fields. The assigned priorities of these 
candidates were later used by the AAOmega fiber-allocation 
software to prioritize between them in cases of conflict. The se¬ 
lected candidates span a range of 0.1-17 days in predicted orbital 
periods and of 0.3-87 mmag in photometric amplitudes. As we 
show below, the confirmed candidates span somewhat narrower 
ranges of 0.3-10days in orbital period and of 0.5-87mmag in 
photometric amplitudes. 

3. AAOmega follow-up observations 

We performed RV follow-up o bservations with the AAOmega 
multi - object spectrograph ( e.g., iLewis et ~aT1 l2002t ISmith et all 
l2004t ISaunders et alJl2004l) at the Anglo-Australian Telescope 
(AAT). The AAOmega two-degree field of view, its three- 
magnitudes dynamic range, and its ability to simultaneously 
record up to 392 spectra are ideal for our purpose. These features 
enabled us to observe the majority of BEER candidates in each 
of the CoRoT fields in a single pointing. We used the AAOmega 


softw are CONFIGURE (e.g., ILewis et all 120021 iMiszalski et al.l 
l2006lPl to optimize the pointing and fibre allocation. Table[l]lists 
for each field the selected pointing ephemeris, the number of 
BEER targets observed, and the total number of science targets 
observed. The observed stars span the range of 12.5-16 in V 
magnitude. 

In total, we observed 281 out of the 481 selected BEER can¬ 
didates. Most of the candidates that were selected for follow-up 
but were not observed are from the LRcOl and LRc02 fields. 
These two fields were observed by CoRoT with two 1.4° x 1.4° 
detectors, as opposed to the LRc03, LRc04, and LRc05 fields 
that were obs erved by CoRoT with only one such detector 
(iMoutou et al.ll2013l) . Nevertheless, to maximize the number of 
priority 1 candidates observed each night, we used only one 
pointing per field. Table [2]lists the coordinates, magnitudes, pho¬ 
tometric ephemeris, and amplitudes of the three BEER effects 
for the candidates observed with AAOmega. For convenience, 
the confirmed binaries are indicated in the rightmost column of 
Table [2] To make the best use of available observing resources, 
we also observed several hundred CoRoT EBs. In this paper, 
however, we report only on observing non-eclipsing BEER can¬ 
didates and on confirming seventy of them. We leave the EBs 
spectra analysis to a separate paper. 


1 Available at http://www.aao.gov.au/science/software/configure 
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Frequency [1/day] Phase 


Fig. 2. BEER light curve analysis of CoRoT 104674562. The panels are the same as in Figure|T] 
Table 2. Coordinates and photometric parameters of the BEER candidates observed at AAOmega" 


CoRoTID 

RA 

(deg) 

Dec 

(deg) 

V 

(mag) 

Orbital 

period 

(day) 

Orbital 

phase 

(HJD-2451545) 

Ellipsoidal 

amplitude 

(ppm) 

Beaming 

amplitude 

(ppm) 

Reflection 

amplitude 

(ppm) 

conf. 

flag 4 

100537909 

290.7709 

1.2476 

13.1 

5.75 

2828.9 

-641 

406 

-216 

1 





0.19 

2.5 

267 

207 

204 


100576007 

290.8286 

1.6783 

15.1 

0.319485 

2833.575 

-902 

227 

-187 

0 





0.000074 

0.017 

27 

21 

44 


100604403 

290.8710 

1.6941 

13.7 

3.094 

2832.39 

-164 

327 

120 

0 





0.021 

0.50 

24 

17 

26 


100637229 

290.9165 

1.4645 

13.8 

8.73 

2829.7 

-675 

480 

-148 

0 





0.17 

1.4 

45 

90 

88 


100677124 

290.9786 

0.9625 

15.7 

0.41720 

2833.533 

-9014 

650 

-1291 

0 





0.00015 

0.026 

67 

581 

420 



(<l) This table is available in its entirety in a machine-readable form in ftp://wise-ftp.tau.ac.il/pub/corotAAO A portion of the table is shown here 
for guidance regarding its form and content. Each line of parameters is followed by a line of uncertainties. 

(4) conf. flag: 0—false alarm; 1—confirmed BEER binary. 


Table 1. Selected pointing ephemeris of the five CoRoT fields, 
and the number of targets observed in each field. 


Field 

Mean RA 
(J2000) 

Mean Dec 
(J2000) 

Science 

targets 

BEER 

targets 

LRcOl 

19:26:25.34 

+01:12:00.5 

211 

41 

LRc02 

18:41:52.44 

+06:37:03.9 

292 

55 

LRc03 

18:32:27.45 

-06:21:52.7 

362 

50 

LRc04 

18:33:50.89 

+08:49:58.9 

359 

68 

LRc05 

18:39:19.29 

+04:28:21.6 

165 

67 


The observations took place on seven consecutive nights, 
starting on August 02, 2012. Between two and five fields 
were observed each night, depending on the available time and 
weather conditions. We acquired six spectra for 231 candidates 
and seven spectra for another 50 candidates in a seven-night 
AAOmega campaign. Tabled lists for each field the heliocentric 
Julian dates (HJDs) of mid-exposure, calculated for the mean 
ephemeris presented in Table [T] 

We used the 1700D grating on the red arm si nce it was re¬ 
ported to give good RV precision (e.g., lLane et ah] 1201 lh . We 
used the 1700B grating on the blue arm since it covers several 
Balmer lines and also enables RV measurements of hot stars. The 
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Table 3. Mid-exposure HJD-2456141 days. 


LRcOl 

LRc02 

LRc03 

LRc04 

LRc05 

1.070956 

— 

0.998684 

— 

— 

2.050900 

1.942608 

1.888283 

1.996221 

2.103633 

3.045416 

2.936347 

2.882764 

2.989653 

3.098236 

4.053797 

3.942036 

3.887762 

3.995172 

4.110621 

5.104335 

4.942736 

4.889214 

4.996094 

5.049920 

6.109585 

5.946059 

5.889688 

6.001038 

6.055340 

— 

6.950716 

6.899660 

7.009869 

7.059933 


nominal spectral coverage is 8342-8842 A on the red arm and 
3803-4489 A on the blue arm, but the actual coverage is smaller 
by up to ~ 60 A and is different for different fibers, depending 
on their position on the detector. 

Observations and data reduction were performed simi larly 
to previously reported works (see e.g. JSebastian et ahll2012l) . In 
short, the observing sequence at each new pointing usually con¬ 
sisted of a flat and two arc frames followed by three or two ex¬ 
posures of 20 or 30 minutes, respectively. 

D ata red uctio n used the dedicated software 2DFDR 
dTavlor et al.1 11996IPI The spectrum from each fiber and 
each subexposure was first normalized by its flat and was 
wavelength calibrated using the arc frames, rebinning the data 
onto the same linear wavelength scale to facilitate combining 
the subexposures. The throughput of each fiber in each subex¬ 
posure was calculated using sky emission lines. An average 
sky spectrum was subtracted from the data using dedicated sky 
fibres. Finally, the subexposures were combined to give a single 
calibrated 60-minute exposure, weighting each subexposure by 
its flux level and rejecting cosmic ray hits. 

Most of the exposures were made under bright sky, so the 
sky subtraction procedure left some residuals of telluric emis¬ 
sion lines in the spectra. In addition, a few pixels in each spec¬ 
trum were affected by bad columns of the detector. We replaced 
the values of the telluric-line residuals and the bad pixels with 
a smoothed version of the same spectrum, which was calculated 
using a moving median filter, 21-pixel wide. Additional 6cr out¬ 
liers were later replaced the same way. 


4. Spectral analysis 

4. 1. Non-composite spectra analysis 

To derive the most precise RVs, we first searched for an optimal 
theoretical template spectrum for each candidate by maximizing 
the cross-correlation values between the candidat e’s obse rved 
spectr a and a set of synthetic Phoenix spectra dHauschildt et al.l 
119991) . calculated on a grid of effective temperature (7^), sur¬ 
face gravity (logg), and metallicity ([m/H]). Line broadening 
due to the projected rotational velocity was added to each syn¬ 
thetic s£ecttum_by_convolvingk wifit a rotational profile G(v) 
(e.g., ISanterne et al.lI 20 T 2 I iGravl [20051) . To account for the in¬ 
strumental broadening of the lines, each synthetic spectrum was 
also convolved with a Gaussian of cr — 20.5 km s _1 . 

Figures [3] and [4] illustrate the data and the optimization pro¬ 
cess. Figure [3] shows co-added AAOmega spectra and cross¬ 
correlation functions (CCFs) of four slowly rotating candi¬ 
dates (rotational broadening < 20 km s _1 ), which have differ¬ 
ent T e ff values and non-composite spectra. The narrow metal 
lines weaken and wide hydrogen lines become stronger with in¬ 
creasing T e ff . Accordingly, the CCF peaks of cool stars (T t fr < 

2 Available at http://www.aao.gov.au/science/software/2dfdr 


H (Balmer) ( a ) 
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Radial Velocity [km/s] 


Radial Velocity [km/s] 


Fig. 3. Co-added AAOmega spectra (panels (a) and (b)) and 
CCFs (panels (c) and (d)) of four BEER candidates with differ¬ 
ent effective temperatures. For clarity, successive spectra were 
shifted upward by 1 and successive CCFs were shifted upward 
by 0.55. CoRoT IDs and estimated temperatures of the four can¬ 
didates are indicated in panel (c). Blue-arm spectra and their 
CCFs are shown on the left (panels (a) and (c)), while red-arm 
spectra and their CCFs are shown on the right (panels (b) and 
(d)). Wavelengths of Call and hydrogen lines are taken from 
iKramida et al.l d2013b . 



Fig. 4. T e ff optimization plot of CoRoT 310199081. Colored dots 
show the maximum-CCF values, cross-correlating its red-arm 
AAOmega spectra with synthetic Phoenix spectra of different 
r eff . Solid black lines show for each exposure a SPLINE interpo¬ 
lation of its maximum-CCF versus 7/n curve. 


7000 K) are narrower than those of hot stars. Figure [4] shows 
the Tetr optimization plot of the coolest candidate presented in 
Fig. [3] The best 7/n for each exposure was estimated from the 
interpolation of its maximum-CCF versus T e s curve. The best 
global 7’eff for each candidate was taken as the weighted mean of 
the different 7’eff values from the different exposures, taking the 
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Table 4. Template parameters of non-composite spectrum 
candidates". 


CoRoT 

ID 

Teff 

(K) 

l°g£ 

(cgs) 

[m/H] 

(dex) 

Rot. broad, 
(km s -1 ) 

conf. 

flag 

100537909 

6500 

5.0 

0.0 

7 

1 

100576007 

5400 

4.0 

-0.5 

1 

0 

100604403 

7000 

4.5 

0.0 

9 

0 

100637229 

6500 

4.5 

0.0 

16 

0 

100677124 

6000 

4.0 

+0.5 

102 

0 


<a> This table is available in its entirety in machine-readable form at 
ftp://wise-ftp.tau.ac.il/pub/corotAAO A portion of the table is shown 
here for guidance regarding its form and content. 


squares of the maximum-CCF values as weights. Surface grav¬ 
ity, metallicity, and projected rotational velocity were estimated 
for each candidate the same way. For each candidate, the near¬ 
est synthetic spectrum to its estimated spectral parameters was 
chosen as its best template. 

The analysis, including RV derivation and orbital solutions, 
was independently performed both for the red-arm and blue-arm 
spectra. However, the results presented from this point onward 
are based on the red-arm spectra alone, since we generally found 
them to have better S/N and stability than the blue-arm spectra. 
For instance, Fig.[3]shows that the blue-arm CCF peaks are wider 
and lower than the red-arm CCF peaks of the same stars. This 
is partially due to the lower resolution and lower S/N of blue- 
arm spectra. Red-arm spectra thus provided better RV precision 
than blue-arm spectra, even for the hottest stars in the sample 
(T eff > 8000 K). 

After optimizing the templates, we derived RVs and errors 
from each spectrum by calculating the CCF with the best tem¬ 
plate. At this stage, we carefully inspected the CCFs, looking 
for a signature of a secondary star, and identified 26 stars that 
present composite spectra (i.e., spectra containing lines of more 
than one star). Table [4] lists the 255 non-composite spectrum 
candidates observed and the template parameters used to derive 
their RVs. For convenience, the rightmost column identifies the 
confirmed BEER single-line binaries (SBIs), 45 of which were 
identified in that list (as explained in Sect. [5}. The measured RVs 
of all non-composite spectrum candidates are given in machine- 
readable form at ftp://wise-ftp.tau.ac.il/pub/corotAAO 

4.2. Composite spectra analysis 

Composite spe c tra were analyzed with TODCOR 
dZucker & Mazehl 119941: iMazeh & Zuckeii 119941) —the two- 
dimensional correlation algorithm. The templates were 
optimized in several steps. First, we searched for the best 
primary template (i.e., the template for the more luminous 
component in the spectrum) similarly to the search performed 
for non-composite spectra. By inspecting the CCFs and the 
derived RVs, we then identified for each candidate the exposures 
made at times when the separation between the primary and 
secondary sets of lines was relatively large. Using only these 
exposure subsets, we optimized the primary- and secondary- 
template T e ff and rotational broadening, also optimizing the 
flux ratio between the primary and the secondary components 
(TODCOR’s a parameter). Metallicity and log g of the secondary 
template (i.e., the template for the less luminous component 
in the spectrum) for each candidate were fixed assuming the 
secondary is a dwarf (log g = 4.5) with the same metallicity as 
its primary. 



8350 8400 8450 8500 8550 8600 8650 8700 8750 8800 

Wavelength [Angstrom] 




-1000 -500 0 500 1000 


Radial Velocity [km/sec] 


Fig. 5. Upper panel: red-arm AAOmega spectrum of the double¬ 
line BEER binary 105962436 from August 05, 2012. Middle: 
TODCOR two-dimensional correlation function for that spectrum, 
using the templates whose parameters are specified in Table 
0 Colored solid lines connect points of equal correlation. The 
black dashed lines run through the two-dimensional correlation 
peak parallel to the primary and secondary RV axes. Bottom 
panel: primary and secondary cuts through the two-dimensional 
correlation peak. 


After optimizing the templates, we used TODCOR to derive 
the primary and secondary RVs and errors from each spectrum, 
fixing a to its best value. As an illustration, Fig.[5]shows a spec¬ 
trum and the TODCOR plots for one of the AAOmega exposures 
of CoRoT 105962436—probably an A5V-G0V double-line bi¬ 
nary (SB2) with a ~ 0.25. The splitting of the Call lines can be 
seen in the upper panel. The middle panel shows the correspond¬ 
ing TODCOR two-dimensional correlation function. The lower 
panel of the figure shows the primary and secondary cuts through 
the two-dimensional correlation function, which run through the 
two-dimensional correlation peak parallel to the primary and 
secondary RV axes dZucker & Mazehll 19941) . The correlation in 
the secondary cut drops only by ~ 0.03 when moving away from 
the peak because we changed the velocity of the secondary tem¬ 
plate, which only contributes ~ 25% of the light, and the primary 
velocity was kept at its best value. Nevertheless, the secondary 
peak is prominent, which means that the RV of the secondary 
star can be measured despite the broad hydrogen lines of the pri¬ 
mary and the relatively low flux ratio. 

Table 0 lists the template parameters and a values used 
for the composite-spectrum BEER candidates. We only use the 
nomenclature ’A’ and ’B’ (both here and later in Table [8]) to 
denote the more and less luminous components, respectively. 
The measured RVs and errors of both components are given in 
machine-readable form at ftp://wise-ftp.tau.ac.il/pub/corotAAO 
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Table 5. Template parameters of the composite-spectrum 
candidates". 


CoRoT 

ID 

Tcfr 

(K) 

l°g£ 

(cgs) 

[m/H] 

(dex) 

Rot. broad. 
(kms~') 

Flux 
ratio (a) 

100688131 A 

6100 

4.5 

0.0 

13 

— 

100688131B 

5900 

4.5 

0.0 

6 

0.62 

100906796 A 

7000 

4.5 

0.0 

60 

— 

100906796B 

5100 

4.5 

0.0 

71 

0.15 

100976101 A 

7400 

4.5 

0.0 

1 

— 

100976101B 

4800 

4.5 

0.0 

16 

0.09 

101177998 A 

7200 

4.5 

0.0 

41 

— 

101177998B 

5000 

4.5 

0.0 

14 

0.17 

103833966 A 

6200 

5.0 

+0.5 

1 

— 

103833966B 

5600 

4.5 

+0.5 

21 

0.27 

(al This table is 

available 

in its entirety in 

machine-readable form at 


ftp://wise-ftp.tau.ac.il/pub/corotAAO A portion of the table is shown 
here for guidance regarding its form and content. 


5. Orbital solutions 

5.1. Confirming BEER binaries with non-composite spectra 

To separate true BEER binaries from false detections (which we 
call false alarms for simplicity, or FAs), and to derive the orbital 
parameters of the true BEER SB Is, we fitted the derived RVs 
of non-composite spectrum candidates with a circular Keplerian 
model. We calculated two x 2 statistics— x^uu an( J f° r the 
null hypothesis (constant RV star) and for the circular orbital so¬ 
lution, respectively. We took the BEER period and phase as pri¬ 
ors by treating them as additional measurements, meaning that 
their squares of residuals, scaled by their error estimates, were 
added to^. in the search for the best fit. For the best-fit orbital 
parameters we also calculated the following F-statistic: 


F = 


(Tnull ~*orb) 


DOF^ - DOF orb 


I A orb 

DOForb ’ 


(7) 


where DOF nu n and DOF orb are the numbers of degrees of free¬ 
dom of the null hypothesis and the orbital solution, respectively. 
Since the only free parameter of the null hypothesis is an RV off¬ 
set, DOF nu n = Nrv - 1, where Nrv is the number of RV points. 
For the orbital solution DOF orb = Nrv + 2 - p olb , where p orb is 
the number of free parameters of the orbital solution, since we 
treated the BEER period and phase as additional measurements. 
For a circular orbit p orb = 4. 

To obtain a good distinction between true BEER SB Is and 
FAs, we used both the ^ n -test and the F-test. The ^ u]1 -test was 
used first to screen all candidates that showed no significant RV 
variability within the observing run. Then the F-test was applied 
to candidates that passed the ^ u|| -test to check the compatibil¬ 
ity of their RVs with a circular Keplerian model at the BEER- 
predicted period and phase. Only candidates that passed both 
tests were considered as confirmed BEER SB Is. 

After visual inspection of the results, particularly of a few 
borderline cases, we chose the critical p-values to be 10 6 and 
0.003 for the X 2 m \\~ an( 3 F-tests, respectively. As a consequence, 
54 out of the 255 non-composite spectrum candidates passed the 
x\ 11 -test, and 45 of them also passed the F-test and were classi¬ 
fied as confirmed BEER SB Is. The nine candidates that passed 
the xl u i \-test but not the F-test might be true variables at a differ¬ 
ent orbital period, or their spectra suffer from some systematics 
causing RV outliers. 

Another possible reason for a true BEER binary to fail our F- 
test is an eccentric orbit. Therefore, we also fitted each RV curve 



Fig. 6. Test statistics p-value histograms of the 255 non¬ 
composite spectrum BEER candidates, solving for their red-arm 
AAOmega RVs. The solid black lines represent a flat distribution 
between 0 and 1 scaled so that its integral is equal the number 
of FAs (210). The dashed black lines mark the critical values 
that were chosen to separate possible BEER SB Is (to the left of 
the lines) from FAs (to the right of the lines). The histograms 
were truncated at 10 7 , so that the leftmost columns contain all 
candidates with /;-values < 10 7 . 


with an eccentric Keplerian model, for which p 0lb = 6. The fact 
we have only 6-7 RV points for most of our candidates makes 
the DOF orb of an eccentric solution as small as 2-3. Since an 
F-test fails at such a low number of DOF, we required a /?-value 
improvement of at least a factor of 10 to prefer the eccentric so¬ 
lution over the circular one. None of the candidates fulfilled this 
requirement, meaning we could not find significant eccentricity 
in any of the confirmed SB Is. 

Figure [6] shows (on a log-scale) /;-value histograms of the 
xl M ~ and F-statistics for the 255 non-composite spectrum BEER 
candidates. The solid black lines represent the expected /;-value 
distributions for RV measurements of constant-RV stars nor¬ 
mally distributed for each star around its RV. The expected and 
observed histograms agree fairly well, particularly for the XnuiC 
statistics. The dashed black lines mark the critical values that 
were chosen to separate possible BEER SB Is from FAs. Only 
candidates found to the left of the lines in both plots were con¬ 
sidered as confirmed BEER SB Is. 

Table [6] lists the orbital parameters of the confirmed BEER 
SB Is. Figure [T2l shows their measured AAOmega RVs and the 
best-fit circular Keplerian model. Their orbital periods span 
a range of 0.4-10days, and their RV semi-amplitudes span a 
range of 6-115 km s~'. Two of the confirmed BEER SBls with 
the smallest RV semi-amplitudes are CoRoT 105659320 and 
101044188—possibly two BDs on a ~ 1 -day period orbit 
around F-G stars. 

5.2. Confirming BEER binaries with composite spectra 

There are three possible scenarios for a composite-spectrum can¬ 
didate. (1) The two components in the candidate’s spectra be¬ 
long to the primary and secondary stars in a short-period SB2 
at the BEER-predicted period and phase. (2) One of the com¬ 
ponents belongs to a binary at the BEER-predicted period and 
phase, while the other component belongs to another star, either 
bound or unbound to the binary (i.e., a diluted BEER binary). 
(3) Neither of the components belongs to the BEER-predicted 
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Table 6. Orbital parameters of the 45 confirmed BEER SB 1 s' 1 . 


CoRoT ID 

P 

(day) 

To 

(HJD-2456141) 

K 

(kms -1 ) 

7 

(kms -1 ) 

F-test 
p -value 

Xtrb 

103924393 

6.851 

5.906 

49.93 

42.90 

1.3E-08 

1.6 


0.021 

0.021 

0.75 

0.60 



105767195 

4.0526 

3.139 

27.99 

2.01 

3.4E-08 

0.7 


0.0022 

0.019 

0.85 

0.60 



105810223 

0.5929447 

3.87271 

115.0 

4.31 

4.0E - 08 

4.5 


0.0000029 

0.00097 

1.4 

0.93 



104674562 

4.5804 

5.895 

30.55 

17.73 

4.5E-08 

0.7 


0.0071 

0.021 

0.96 

0.66 



310222016 

2.08218 

5.725 

29.3 

-30.61 

4.9E - 08 

3.8 


0.00047 

0.020 

1.7 

0.65 



310210722 

6.975 

7.165 

47.7 

11.02 

1. IE - 07 

7.6 


0.020 

0.025 

1.2 

0.78 



105818861 

1.301182 

6.0599 

57.5 

9.82 

1.2E-07 

3.5 


0.000054 

0.0037 

1.0 

0.73 



101029997 

3.57665 

3.361 

41.42 

36.17 

1.3E-07 

2.7 


0.00044 

0.011 

0.88 

0.59 



105844488 

0.554937 

2.9081 

90.2 

-7.7 

1.4E-07 

2.4 


0.000095 

0.0021 

2.1 

1.5 



310215106 

3.23943 

4.099 

33.58 

-0.07 

1.7E-07 

7.4 


0.00032 

0.016 

0.91 

0.69 



310176634 

1.26221 

5.877 

17.1 

-26.82 

3.5E-07 

0.5 


0.00058 

0.013 

1.1 

0.81 



105336757 

2.356960 

6.059 

32.78 

-7.64 

3.8E-07 

3.0 


0.000096 

0.011 

0.86 

0.64 



105712106 

3.4521 

7.165 

47.04 

-2.58 

7.7E - 07 

10.4 


0.0020 

0.010 

0.78 

0.59 



105403147 

2.9136 

5.876 

35.37 

22.14 

1. IE - 06 

7.0 


0.0019 

0.012 

0.79 

0.59 



105706604 

2.26769 

6.978 

32.4 

-17.00 

2.2E - 06 

4.8 


0.00017 

0.012 

1.2 

0.69 



105714214 

3.1309 

1.950 

36.57 

-10.49 

2.4E - 06 

10.0 


0.0051 

0.012 

0.82 

0.60 



104295292 

2.8685 

2.751 

26.09 

7.83 

2.8E - 06 

5.1 


0.0043 

0.015 

0.86 

0.60 



110567660 

1.238837 

3.9782 

49.5 

-8.1 

2.9E - 06 

6.0 


0.000091 

0.0054 

1.5 

1.0 



101177265 

0.6419480 

4.0489 

36.4 

35.2 

3.0E - 06 

3.4 


0.0000050 

0.0038 

1.5 

1.0 



101014035 

0.4596816 

4.0441 

55.3 

1.9 

4. IE - 06 

7.6 


0.0000015 

0.0021 

1.7 

1.2 



100889978 

10.160 

0.734 

27.38 

7.84 

5.4E-06 

9.4 


0.026 

0.093 

0.80 

0.98 



105661774 

2.22043 

2.094 

31.3 

26.2 

6.2E - 06 

2.0 


0.00013 

0.023 

2.3 

1.4 



105618890 

2.6219 

1.774 

39.5 

21.24 

6.5E - 06 

11.0 


0.0012 

0.013 

1.1 

0.81 




Note: each line of parameters is followed by a line of uncertainties. 


binary. We classified a composite-spectrum BEER candidate as 
confirmed if the RVs of at least one of its components was com¬ 
patible with a Keplerian model at the BEER-predicted period 
and phase. 

To assign the correct scenario to each composite-spectrum 
candidate and (in case of a true BEER binary) to derive its orbital 
parameters, we separately fitted a circular Keplerian model to its 
primary and secondary RVs and also calculated the /t-value of 
the F-test in Equation[7| Candidates in which both component’s 
SBl-model got ap-value < 0.001 were considered as confirmed 
BEER SB2s. If the SB1 model of only one of the components 
got a p-value < 0.001, the candidate was considered as a diluted 
BEER binary. If both component’s SBl-model got a p-value > 
0.001, we considered the candidate as an FA. Similarly to non¬ 
composite spectrum candidates, the selected critical p-value of 


0.001 originates in a visual inspection of the results, particularly 
of a few borderline cases. 

In our sample of 26 composite-spectrum candidates we 
found 18 BEER SB2s, 7 diluted BEER binaries, and one FA 
(CoRoT 310186704). For the confirmed BEER SB2s we then 
also fitted circular SB2 Keplerian model for the two sets of RVs 
together. Together with the 45 confirmed BEER SB Is, we have 
thus confirmed 70 new non-eclipsing BEER binaries. 

In addition to a circular model, we also fitted each RV curve 
with an eccentric Keplerian model. Requiring an improvement 
of a factor 10 in the F-test />value, we found no SB2s that 
show measurable eccentricity. We found two diluted-binary can¬ 
didates (104626523 A and 103833966 A) to have slightly eccen¬ 
tric orbits. However, the low eccentricities (~ 0.3) found for 
these two cases might also be spurious or at least inflated (e.g., 
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Table 6. Continued. 


CoRoTID 

P 

(day) 

T 0 

(HJD-2456141) 

K 

(kms -1 ) 

7 

(kms -1 ) 

E-test 
p -value 

Xtrb 

105802223 

1.414011 

3.826 

18.0 

34.40 

7.6E - 06 

2.0 


0.000040 

0.014 

1.2 

0.85 



100880613 

0.893171 

1.9956 

36.0 

59.9 

8.0E - 06 

8.3 


0.000018 

0.0069 

1.3 

1.1 



105760939 

1.53138 

3.847 

22.7 

36.56 

8.3E - 06 

2.6 


0.00042 

0.012 

1.4 

0.89 



104113878 

2.22939 

3.085 

15.70 

-25.11 

2.5E - 05 

6.5 


0.00024 

0.024 

0.97 

0.59 



103782315 

5.372 

6.737 

15.9 

-46.68 

3.IE - 05 

2.9 


0.012 

0.070 

1.3 

0.90 



310169750 

0.6321732 

2.8879 

40.5 

-3.8 

3.4E - 05 

8.0 


0.0000065 

0.0058 

2.0 

1.5 



103922738 

7.901 

7.92 

6.2 

12.64 

5.3E-05 

0.5 


0.037 

0.19 

1.4 

0.87 



104648865 

1.446104 

5.045 

32.6 

-21.2 

6.3E - 05 

0.8 


0.000097 

0.041 

5.9 

4.3 



104536524 

3.12477 

2.108 

10.8 

60.34 

6.7E - 05 

2.6 


0.00100 

0.052 

1.1 

0.79 



101058035 

0.7647344 

6.1885 

21.3 

10.57 

8.6E - 05 

9.1 


0.0000057 

0.0079 

1.3 

0.95 



105378453 

3.2362 

5.099 

19.77 

-14.99 

1.5E-04 

23.1 


0.0075 

0.021 

0.83 

0.58 



105597526 

0.920035 

2.934 

39.0 

-55.7 

2.0E - 04 

7.1 


0.000029 

0.023 

3.7 

5.8 



104667709 

0.6984611 

3.0336 

25.1 

-34.8 

3.5E - 04 

13.1 


0.0000056 

0.0071 

1.9 

1.1 



105154613 

4.587 

1.459 

9.36 

38.69 

5.0E - 04 

8.6 


0.026 

0.064 

0.88 

0.59 



100537909 

4.8190 

2.534 

24.41 

31.81 

5.6E - 04 

67.6 


0.0094 

0.023 

0.83 

0.56 



105659320 

0.706141 

4.986 

6.91 

-1.18 

6.7E - 04 

4.7 


0.000026 

0.018 

0.94 

0.72 



101044188 

1.368287 

4.969 

8.19 

39.32 

7. IE - 04 

6.9 


0.000075 

0.021 

0.94 

0.62 



104279119 

10.161 

7.63 

14.8 

12.5 

7.3E - 04 

15.2 


0.051 

0.24 

1.5 

1.7 



104598628 

2.75207 

4.264 

19.1 

96.2 

1.3E-03 

6.5 


0.00061 

0.056 

2.8 

1.8 



105164611 

3.50481 

4.532 

23.8 

4.0 

1.4E-03 

9.0 


0.00066 

0.056 

2.9 

1.8 



100851348 

0.854818 

4.027 

13.5 

19.0 

1.7E-03 

4.0 


0.000012 

0.023 

2.5 

2.0 



105472536 

0.4066543 

6.8778 

34.0 

-36.8 

2.IE - 03 

4.8 


0.0000092 

0.0050 

6.4 

2.4 




iLucv & Sweenevil 19711) . The preference for an eccentric orbital 
solution might also be a result of systematic RV errors caused by 
the presence of the second component in the spectra. 

Table [7] lists the orbital parameters of the confirmed BEER 
SB2s, sorted by ascending E-test p-value of the primary. Table 
[8] lists the orbital parameters of the variable components in the 
diluted BEER binaries, sorted by ascending E-test p-value. 

Figure [T3l shows the measured AAOmega RVs and the best- 
fit Keplerian models of the confirmed BEER SB2s listed in Table 
|7] Figure[14]shows the measured AAOmega RVs and the best-fit 
Keplerian models of the variable components in the confirmed 
BEER diluted binaries listed in Table [3 

Figure [15] shows the phase-folded and binned light curves of 
all 70 confirmed BEER binaries, together with the best-fit circu¬ 
lar BEER model. For convenience, the order of the plots in Fig. 
n~5lis the same as in Figs. [T2l4T4l 


6. Performance of the BEER search algorithm 

We are now in a position to evaluate the performance of the 
BEER algorithm in detecting short-period binaries in the light 
curves of CoRoT long runs. This is possible in view of the large 
sample of confirmations and FAs that are reported here. 


6.1. BEER-model priority, M 2 sin i, and period 

We start by considering the priority classes that we manually 
assigned to candidates during the visual inspection stage. The 
left-hand side of Table [3 lists the number of binary confirma¬ 
tions and FAs per priority class within our sample. As expected, 
BEER did not perform that well with its priority 3 candidates. 
They were not considered to be good candidates in the first place 
and were included in the observational campaign only due to the 
availability of fibers on the AAOmega spectrograph. Therefore, 
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Table 7. Orbital parameters of the 18 confirmed BEER SB2s". 


CoRoTID 

P 

(day) 

To 

(HJD-2456141) 

Ki 

(kms -1 ) 

k 2 

(kms -1 ) 

7 

(km s -1 ) 

7-test 
p-value A 

7-test 
p-value B 

Xorb 

105962436 

1.8020 

5.917 

83.5 

133.8 

27.74 

5.6E-09 

8.0E - 05 

20.3 


0.0062 

0.012 

1.2 

3.0 

0.78 




106024478 

3.013 

4.1436 

65.25 

66.2 

-22.33 

2.0E - 07 

5.4E - 08 

9.9 


0.012 

0.0068 

0.98 

1.1 

0.57 




105649738 

3.604 

2.7189 

87.6 

87.7 

-7.47 

2.0E - 07 

2.2E - 07 

8.4 


0.010 

0.0064 

1.1 

1.2 

0.57 




105928477 

1.6509 

4.9285 

74.3 

74.7 

-2.18 

2.3E - 07 

1.6E-06 

3.1 


0.0048 

0.0063 

1.2 

2.3 

0.80 




100688131 

6.984 

2.518 

38.12 

43.20 

-13.84 

1.7E-06 

2.3E - 06 

25.1 


0.098 

0.023 

0.77 

0.95 

0.47 




310198235 

1.713 

6.873 

99.7 

113.4 

7.55 

2.0E - 06 

6.3E - 07 

74.8 


0.011 

0.019 

2.8 

3.7 

0.61 




100976101 

4.738 

0.856 

67.22 

113.0 

-9.00 

2.3E - 06 

2.2E - 05 

68.7 


0.025 

0.018 

0.73 

2.9 

0.55 




104369937 

6.093 

5.699 

56.37 

72.8 

-16.58 

2.6E - 06 

2.2E - 04 

54.7 


0.061 

0.020 

0.81 

1.7 

0.54 




105963904 

6.048 

5.062 

67.26 

81.7 

-36.46 

2.8E - 06 

4.7E - 05 

99.9 


0.044 

0.010 

0.76 

1.5 

0.53 




104181232 

4.138 

4.025 

58.15 

83.2 

-0.64 

4.0E - 06 

1.9E-04 

34.1 


0.027 

0.012 

0.87 

3.5 

0.65 




105506915 

0.6542 

3.083 

25.50 

123.4 

-72.51 

2.0E - 05 

4.8E - 05 

22.8 


0.0048 

0.016 

0.98 

5.0 

0.96 




310136399 

6.146 

2.826 

73.2 

76.0 

-7.10 

3.3E-05 

2.0E - 05 

96.8 


0.098 

0.027 

1.9 

1.7 

0.97 




101177998 

1.6103 

4.0305 

26.5 

156.8 

40.59 

2.8E - 04 

2.4E - 04 

32.7 


0.0044 

0.0056 

1.7 

4.1 

0.91 




310212616 

0.35918 

3.8372 

96.7 

149.0 

-12.0 

3.4E-04 

1.1E-05 

105.1 


0.00047 

0.0016 

3.6 

3.7 

2.0 




103838038 

6.238 

4.758 

42.98 

44.0 

-8.99 

3.8E - 04 

4.4E - 05 

192.2 


0.061 

0.016 

0.88 

1.0 

0.48 




310173237 

0.29882 

2.0277 

85.5 

132.2 

-38.6 

4.7E - 04 

3.9E - 04 

177.7 


0.00033 

0.0025 

3.9 

4.6 

2.1 




104432741 

0.43955 

4.0074 

86.5 

166.1 

9.7 

6.0E - 04 

1.3E-05 

29.8 


0.00053 

0.0018 

2.7 

4.1 

2.0 




105583867 

0.8567 

2.049 

54.4 

141.5 

10.2 

6. IE - 04 

1.4E-04 

47.2 


0.0037 

0.017 

2.7 

7.0 

2.0 





(a> Each line of parameters is followed by a line of uncertainties. 


Table 8. Orbital parameters of the variable component in the seven confirmed diluted BEER binaries". 


CoRoTID 

P 

(day) 

j'b 

(HID-2456141) 

K 

(kms -1 ) 

7 

(kms -1 ) 

7-test 
p -value 

Xtrb 

310205770B 

5.543 

5.218 

56.8 

3.51 

4.9E - 07 

12.6 


0.014 

0.019 

1.5 

0.94 



310193013 B 

0.2801751 

5.8446 

193.3 

-0.3 

1.5E-05 

32.0 


0.0000021 

0.0013 

3.9 

3.3 



104791410B 

1.04457 

6.757 

108 

0 

2.0E - 05 

5.2 


0.00051 

0.021 

18 

18 



104626523 A 

7.354 

2.99 

28.2 

4.7 

LIE - 03 

0.1 


0.052 

0.22 

1.6 

1.5 



105423352B 

4.712 

3.668 

56.5 

-24.0 

2.6E - 04 

10.8 


0.018 

0.040 

2.5 

1.8 



100906796B 

0.8974124 

5.0351 

180.0 

-6.9 

3.5E - 04 

39.7 


0.0000046 

0.0043 

4.8 

3.2 



103833966 A 

10.80 

7.75 

72 

1 

4.0E - 03 

0.4 


0.11 

0.92 

47 

22 




Each line of parameters is followed by a line of uncertainties. 

For the eccentric binaries 104626523 A and 103833966 A the third column shows the periastron time. The eccentricity of 104626523 A is 
0.33 ± 0.09 and its longitude of periastron is 242 ± 15 deg. The eccentricity of 103833966 A is 0.32 ± 0.16 and its longitude of periastron is 
121 ±40 deg. 


we decided to ignore priority 3 candidates in our performance Two other parameters that determine the success of the 
analysis, and only concentrated on priority 1 and 2 candidates. BEER detection are the secondary mass and the orbital period 
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Table 9. Confirmations and FAs per priority class by applying different filters to the sample 11 . 


Filter: 

All candidates 

Mi sin i > 0.25A/ 0 

Priority 

Confirmations 

FAs 

Confirmations 

FAs 

1 

30 

26 

30 

18 

2 

35 

108 

35 

71 

3 

5 

77 

5 

70 


Mi sin i > 0.25 M e and spectral 
type earlier than G7 for P > 1.4 day 
Confirmations FAs 


30 

34 

4 


17 

54 

35 


(a) See text for the justifications of the chosen filters. 
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Fig. 7. BEER-model M 2 sin i as a function of the BEER pe¬ 
riod. Large and small circles represent priority 1 and 2 tar¬ 
gets, respectively. Green circles represent confirmations, while 
red circles represent FAs. Open circles represent candidates of 
spectral type later than G7. The horizontal dashed line marks 
M 2 sin i = 0.25 M 0 , below which there are no confirmed tar¬ 
gets. The vertical dashed line marks P - 1.4 days, to the right 
of which there are 18 FAs of spectral type later than G7, but 
only one confirmed binary (see Sect. 16.2k 


because t he BEER amplitudes m ostly depend on these two pa¬ 
rameters dFaigler & Mazehl201 lh . We therefore plot in Fig.[7]the 
BEER-model M 2 sin i as a function of the photometric period for 
priority 1 and 2 targets. As expected, shorter-period modulations 
with larger BEER-model M 2 sin i have higher chances of being 
true BEER modulations than FAs. For instance, there are no con¬ 
firmations with BEER-model M 2 sin i < 0.25M Q . The vertical 
dashed line in Fig. [7] marks P = 1.4 days, and it is explained 
next. 

From this stage onward we continue the performance anal¬ 
ysis while ignoring all targets with BEER-model M 2 sin i < 
0.25 M a . The summary of the remaining targets as a function 
of their priority is listed in the central columns of Table [9] 

6.2. Spectral type 

Another parameter that can influence the success of the BEER 
detection is the candidate’s spectral type. To check this option, 
we plot in Fig. [8] the BEER score as a function of the spectral 
type for priority 1 and 2 candidates with BEER-model M 2 sin i > 
0.25M q . The spectral type is taken from EXODAI0, and it 
was ob tained using the SED analysis described bv lDeleuil et al.l 
(120091) . The advantage of using EXODAT is its availability re- 

3 cesam.oamp.fr/exodat 


Fig. 8. BEER score as a function of the EXODAT spectral type 
for candidates with M 2 sin i > 0.25M Q . The symbols are the 
same as in Fig. [7] The vertical dashed line marks a spectral type 
of G7. 


gardless of any follow-up observations. Fig. [8] shows that there 
is a larger fraction of FAs for late-type candidates than for early- 
type candidates. Particularly if we draw a line at a spectral type 
of G7 (T e ff ~ 5600 K), we find that the fraction of FAs is ~ 54% 
(69/128) to the left of that line and ~ 77% (20/26) to the right 
of that line. 

In Fig. |7] we also mark with open circles candidates of spec¬ 
tral type later than G7. A close examination of these cool can¬ 
didates, with BEER-model M 2 sin i > 0.25M 0 , reveals that the 
confirmed binaries differ from the FAs in yet another way. By 
drawing a line at P — 1.4 days, we find most of the cool FAs 
(18/20) to the right of that line, and most of the confirmed cool 
binaries (5/6) to the left of that line. 

To explain this phenomenon, we plot in Fig. [3 the spectral 
parameters of the observed candidates, which were derived in 
Sect. [4] Since spectral type is a proxy of T e g, we place //.fi¬ 
at the abscissa. Similarly to Fig. [3 there is an excess of cool 
FAs (T e ff < 5600 K), with just a few confirmed binaries in that 
temperature regime. In addition, cool FAs seem to constitute a 
distinct sample of slowly rotating stars with lower gravity and 
lower metallicity. Even though our spectral-paramete r measu re- 
ment technique is prone to systematic biases fe.«.. iTorres et al.l 
12012h . such a strong bimodality suggests that most cool FAs are 
possibly red giant stars. 

We propose that the main reason that red giants introduce 
false candidates with photometric periods of >1.4 da y is related 
to solar-like oscillations. Using Kepler data. lMosser et al.l ( 2013 ) 
have shown that solar-like oscillations of 1-2 M G red giants have 
frequencies of 1-10//Hz (periods of 1.2-12 days), amplitudes 
(A max ) of 0.1-1 mmag, and that the amplitude increases with the 
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Teff [K] 


Fig. 9. Spectral parameters derived in Sect. [4] Confirmed BEER 
binaries are represented in green, while FAs are represented in 
red. For SB2s, only the spectral parameters of the primary are 
shown. For diluted BEER binaries, only the spectral parameters 
of the variable component are shown. The upper panel shows the 
two :r eff histograms. The lower three panels show the 7/n-log g, 
Teff-metallicity, and /' c n-rotational broadening scatter plots. The 
vertical dashed lines mark T e g — 5600 K—the approximate 
temperature of a G7V star. The horizontal dashed line in the 
lower panel indicates a rotational broadening of 20kms _1 , be¬ 
low which the measured values are unreliable as a result of the 
medium resolution of the spectrograph (see Sect.[4]i. 


period as a power law. For BEER bi naries, the pho tome t ric am- 
plitude decreases with the period (iZucker, Mazeh. & Alexander! 
l2007t) . However, it appears that for periods of about l-10days, 
these two phenomena might have similar photometric ampli¬ 
tudes. A large fraction of the stars observed by CoRoT are 
indeed giants, par ti cularly in the CoRo T-center fields (e.g., 
iDeleuil et al.l l2009t iGazzano et al.l l2010h . Although solar-like 
oscillations are semi-regular in nature, given the typical length 
and S/N of long-run CoRoT data, BEER might have interpreted 
the variability of some red giants as an indication for the pres¬ 
ence of a short-period companion. 

To check whether filtering out red candidates with periods of 
>1.4 days improves BEER’s performance, we list in the right- 
hand side of Table [9] the number of binary confirmations and 
FAs per priority class, after removing candidates with BEER- 
model Mi sin i < 0.25M Q , and also candidates of spectral type 
later than G7 with photometric periods of > 1.4 days. These two 
filters applied to the candidate list lowers the fraction of FAs 
to ~ 1/3 for priority 1 candidates, and to ~ 2/3 for priority 2 
candidates. 




Fig. 10. Upper panel: False-alarm probability as a function of 
the BEER score threshold. Lower panel: The detection proba¬ 
bility as a function of the BEER score threshold. The dashed 
lines show the probabilities for the whole sample of 281 candi¬ 
dates observed, while the solid lines show the probabilities for 
the remaining 135 candidates, after applying the three filters dis¬ 
cussed in the text (see also the first two rows of the right-hand 
side of Table [9]). The plot was truncated at a score threshold of 
0.85 since only 8 candidates got a higher score. 


6.3. BEER-model score 


To evaluate the BEER performance as a function of its score, 
we counted the targets with scores higher than some threshold T 
and obtained the detection probability Pq(T ) and the false-alarm 
probability Pf A {T) of the algorithm as 


Pd = 


Number of confirmations with score > T 
T otal number of confirmations 


( 8 ) 


p Number of false alarms with score > T 

Total number of targets with score > T 

To illustrate the algorithm performance, Fig.[T0]plots Pq and 
Ppa as a function of T for two sets of targets—the whole sample 
of 281 candidates observed, and the remaining 135 candidates, 
after filtering out priority 3 candidates, candidates with BEER- 
model Mi sin i < 0.25M e , and candidates of spectral type later 
than G7 with photometric periods of > 1.4 day. The addition 
of these three filters lowers P fa and raises If , for any given T, 
hence improves the algorithm performance. 

Estimating the two probabilities for any selected threshold 
level might be useful in estimating the results of future CoRoT- 
based BEER searches. If such a search would yield N targets 
with scores higher than some predefined threshold T, then using 
the Pq and P FA that correspond to T, we expect to have 

N c = N(\ - P FA ) , (10) 


Nfa = NPfa , and 


( 11 ) 


N C c = m - PFA) 
Pd 


( 12 ) 


12 















































Tal-Or, Faigler, and Mazeh: Seventy new CoRoT BEER binaries 


Table 10. Spectral parameters, primary mass, and the assigned Table 11. Spectral parameters, primary mass, and mass ratio of 
mass ratio of the 45 confirmed BEER SB Is. the 18 confirmed BEER SB2s. 


CoRoT 

ID 

7eff 

(K) 

l°g£ 

(cgs) 

[m/H] 

(dex) 

Primary mass 
(Mo) 

Mass ratio 
(M 2 /M,) 

103924393 

6500 

4.6 

-0.1 

1.20 + 0.16 

0-69/op 

105767195 

6400 

4.6 

-0.4 

1.07 + 0.13 

0-28/°g 

105810223 

6500 

3.8 

0.0 

1.59 + 0.38 

0-62*|j} 

104674562 

6500 

4.6 

-0.2 

1.15 + 0.15 

0-32!| 

310222016 

6800 

4.5 

-0.1 

1.28 + 0.17 

0.21*}}® 

310210722 

7800 

3.9 

0.1 

2.02 + 0.46 

0.52^3? 

105818861 

6400 

4.6 

-0.2 

1.12 + 0.15 

0 4i781 

0 06 

101029997 

6600 

4.6 

-0.1 

1.20 + 0.15 

0 41*°- 35 
w 0 06 

105844488 

5000 

5.5 

+0.3 

1.01 +0.23 

0.53*°-®° 

310215106 

6300 

3.9 

-0.2 

1.40 + 0.32 

0.28*} 

310176634 

5000 

5.5 

+0.5 

1.07 + 0.26 

o.io! $ 

105336757 

6500 

4.6 

-0.2 

1.16 + 0.15 

°- 26+ 004 

105712106 

6300 

4.6 

-0.0 

1.15 + 0.15 

0.48*}}} 

105403147 

6700 

4.5 

-0.2 

1.24 + 0.17 

0.30*}}? 

105706604 

6300 

4.5 

-0.1 

1.15 + 0.16 

0-26+S 

105714214 

6200 

5.1 

0.2 

1.19 + 0.17 

0.33*}}? 

104295292 

6200 

4.7 

-0.0 

1.12 + 0.14 

0 22* 
v7.z-z._oo3 

110567660 

6500 

4.5 

0.1 

1.24 + 0.17 

0.32*® 

101177265 

6700 

4.2 

0.0 

1.42 + 0.27 

0.17*} 

101014035 

6100 

4.2 

-0.0 

1.24 + 0.24 

0.25*}}} 

100889978 

6300 

5.0 

-0.2 

1.11 +0.14 

0.39*}}} 

105661774 

8900 

4.0 

-0.5 

1.96 + 0.40 

0.20*}}} 

105618890 

6300 

4.4 

0.0 

1.20 + 0.18 

0.34*}| 

105802223 

6400 

4.6 

0.0 

1.20 + 0.16 

an*}! 

100880613 

7300 

4.5 

+0.4 

1.60 + 0.22 

0.18*} 

105760939 

6200 

4.2 

0.1 

1.28 + 0.23 

q i4+8.¥§ 

1 -0.02 

104113878 

6700 

4.7 

0.1 

1.27 + 0.15 

0 . 11 *°°/ 

103782315 

5900 

4.6 

+0.3 

1.12 + 0.15 

°- i6+ }m 

310169750 

6500 

3.6 

+0.4 

1.96 + 0.52 

0 17* 8 ® 

1 ' -0.03 

103922738 

6000 

4.7 

0.1 

1.09 + 0.14 

0 07 +0 -° 5 

W - W/ _Q.01 

104648865 

8700 

3.3 

-0.4 

2.89 ± 0.83 

0.154s 

104536524 

6400 

4.4 

0.1 

1.25 + 0.20 

0.08*}°® 

101058035 

6500 

4.7 

0.1 

1.23 + 0.15 

0.11*} 

105378453 

6300 

5.0 

-0.1 

1.13 + 0.14 

0.17*} 

105597526 

6200 

3.0 

-0.0 

2.47 ± 0.78 

0 17* 8 ® 

W,A/ 004 

104667709 

6200 

3.6 

-0.0 

1.70 + 0.46 

0 n*°.° 7 

105154613 

6300 

4.8 

-0.2 

1.08 + 0.13 

0.09*°°® 

100537909 

6500 

5.0 

-0.1 

1.19 + 0.15 

0.24*}}} 

105659320 

6500 

4.6 

-0.1 

1.16 + 0.15 

0 03 + : ^ 

w.uz>_ooi 

101044188 

5900 

4.8 

0.1 

1.07 + 0.13 

0.05*}° 3 

104279119 

6200 

5.1 

+0.3 

1.24 + 0.18 

0.18*} 

104598628 

6500 

3.7 

0.2 

1.79 + 0.46 

0.13*}f 

105164611 

8500 

4.5 

-0.1 

1.75 + 0.23 

0.18! | 

100851348 

6300 

3.9 

+0.3 

1.58 + 0.36 

0.06*°°! 

105472536 

5500 

3.9 

-0.0 

1.20 + 0.28 

0 1 4+8:93 

_-M3_ 


where Nc is the number of expected confirmations within the 
N candidates, Nfa is the number of expected FAs within the N 
candidates, and Ncc is the estimated number of binaries that can 
be discovered by BEER in the original sample of CoRoT targets. 
For instance, using the solid lines in Fig. [TOl for T = 0.6 we get 
Pfa ~ 1/3 and Pd ~ 1/2. This means that, for T = 0.6, we 
expect about two thirds of the selected sample to be true BEER 
binaries and the total number of BEER binaries in the original 
sample of CoRoT targets to be about (4/3 )N. 

7. Mass ratio and orbital period distribution of the 
BEER CoRoT sample 

To discuss the mass ratio and period distribution of the sam¬ 
ple of confirmed BEER binaries, we wish to plot an estimate 


CoRoT 

ID 

Teg 

(K) 

logg 

(cgs) 

[m/H] 

(dex) 

Primary mass 
(Ms) 

Mass ratio 
(M 2 /Mj) 

105962436 

8200 

4.4 

-0.2 

1.68 + 0.26 

0.62 ± 0.04 

106024478 

6000 

3.2 

-0.1 

2.04 ± 0.63 

0.99 ± 0.02 

105649738 

5800 

3.1 

+0.3 

2.26 + 0.71 

1.00 + 0.02 

105928477 

6300 

4.4 

-0.0 

1.19 + 0.18 

0.99 ± 0.07 

100688131 

6100 

4.7 

0.1 

1.14 + 0.14 

0.88 ± 0.02 

310198235 

6500 

3.9 

+0.3 

1.66 + 0.39 

0.88 + 0.12 

100976101 

7400 

4.7 

0.1 

1.48 + 0.16 

0.60 ± 0.04 

104369937 

6400 

5.0 

0.1 

1.18 + 0.14 

0.77 ± 0.03 

105963904 

6100 

4.7 

0.1 

1.11 ±0.14 

0.82 ± 0.02 

104181232 

6500 

5.0 

-0.0 

1.19 + 0.15 

0.70 + 0.10 

105506915 

6500 

3.8 

0.0 

1.62 + 0.40 

0.21+0.04 

310136399 

5900 

3.3 

+0.5 

2.04 ± 0.60 

0.96 ± 0.06 

101177998 

7200 

4.3 

-0.0 

1.50 + 0.24 

0.17 + 0.03 

310212616 

6100 

4.5 

+0.3 

1.22 + 0.17 

0.65 + 0.11 

103838038 

6300 

4.5 

-0.1 

1.13+0.16 

0.98 ± 0.03 

310173237 

5200 

4.5 

0.2 

0.92 + 0.14 

0.65 + 0.15 

104432741 

5800 

2.9 

+0.4 

2.74 ± 0.90 

0.52 ± 0.07 

105583867 

7300 

3.8 

0.0 

1.92 + 0.47 

0.38 + 0.14 


of the mass ratio of each system as a function of its orbital 
period. For SB2s this is straightforward, since the mass ratio 
(q = M 2 /M 1 = K 1 /K 2 ) can be measured directly from the or¬ 
bital parameters. For SB Is, however, the mass ratio depends not 
only on the orbital parameters, but also on the primary mass and 
orbital inclination, and can be found using the relation 

sin 3 i)q 3 -q 1 -2q-} = 0, (13) 

where M\ is the estimated mass of the primary, / is the mass 
function derived analytically from the orbital parameters, and i 
is the inclination. 

To estimate M\ of the_confirmed binaries, we used the em¬ 
piric relations given in [Torres et~aH (12010l) . which express the 
stellar mass and radius in terms of its observed spectral param¬ 
eters. As an input, we used the spectral parameters that were 
derived for the primary component of the confirmed SB Is and 
SB2s in Sect. [4] To estimate the mass uncertainties, we took 
equal uncertainties of 300 K in T e g, 0.4 dex in log g, and 0.3 dex 
in [m/H] to all primary stars. For most of the observed stars 
these uncertainties are larger than the scatter of the best spectral 
parameters between consecutive exposures (see for instance Fig. 
4]>, but taking into account the possible systematic errors (e.g., 
Torres et aJ1l20l 2l) . they are prob ably reasonable. T he intrinsic 
scatter from the empiric relations (ITorres et al.ll2010t) was added 
in quadrature to the mass uncertainties. 

To assign an inclination for each SB1 we considered three 
options—to use the inclination estimated from its light curve 
by BEER, to derive its inclination distribution using the algo¬ 
rithm developed by iMazeh & Goldberg! ( 1992 ). or to assume 
an isotropic inclination distribution. While using the BEER- 
model sin / estimates could sound appealing, we decided not 
to use them because we do not know their actual uncertain¬ 
ties, including possible systematic biases. Deriving incl i nation 
distributions using the algorithm of iMazeh & Goldberg! dl992 ) 
would probably be the correct way to proceed, but it is beyond 
the scope of this paper. We therefore, somewhat arbitrarily, as¬ 
signed a value of sin i = 0.866 1 taking the median of an 
isotropic inclination distribution and the confidence li mits to 
cover the central 68.3% of t he distribution (e.g.. iHo & Turned 
1201 ltlEonez & Jenkinsll2012h . 


13 































Tal-Or, Faigler, and Mazeh: Seventy new CoRoT BEER binaries 



Period [day] 


Fig. 11. Assigned mass ratio of the 45 confirmed SB Is and esti¬ 
mated mass ratio of the 18 confirmed SB2s versus orbital period. 
For the SB Is we assumed an isotropic inclination distribution 
(see text). Bottom panel: period histogram. The_solid red line 
shows the log-normal fit reported by IRaghavan et all (1201 Ol) to 
the period distribution of stellar binaries, scaled to best fit the 
histogram. Right panel: mass-ratio histogram. 


Table [TOllists the spectral parameters, primary mass, and the 
assigned mass ratio of the confirmed BEER SBls. TableflTI lists 
the primary spectral parameters, primary mass, and the estimated 
mass ratio of the confirmed BEER SB2s. For convenience, the 
order of Tables ITOl and ITT! is the same as of Tables [6] and [ 7 ] We 
did not estimate the masses of the seven diluted BEER binaries 
since their spectral parameters might have been biased by the 
presence of the third star. 

Figure QT| shows the assigned mass ratio of the 45 confirmed 
SBls and the estimated mass ratio of the 18 confirmed SB2s as 
a function of their orbital period. The period and the mass-ratio 
histograms are plotted as well. The solid red line in t he per iod 
histogram shows the log-normal fit of lRaghavan et al.l (1201 Ol) to 
the period distribution of stellar binaries in the solar neighbor¬ 
hood, scaled to best fit the histogram. 

One feature in Figure [TT] is that for periods of > 1 day, there 
seems to be a transition fro m SBls to SB2s at q ~ 0. 6. A similar 
transition was pointed outby lHalbwachs et al.l ( 2003 1 in their un¬ 
biased sample of solar-like spectroscopic binaries. In our sample 
this transition can be explained by the mass-luminosity relation 
of normal stars t hat is expected to be L ~ M 4 (e.g.. lTorres et alJ 
l20lQ[lCoxl20()0L p. 382). Since the lowest significant a we could 
detect was ~ 0.1 (see Table 0, we did not expect detached bina¬ 
ries with q < 0.6 to be classified as SB2s. For periods of < 1 day 
this transition is less evident, since at such short periods some 
of these binaries cou ld experience mass transfer via Roche-lobe 
overflow (e.g JEggletonl 19831) . and the mass-lumin osity relation 
deviates from that of isolated normal stars (e.g.. lBattenlll973L p. 
154). 

The period histogram clearly shows a rise in the number 
of binaries per AlogP from ~ 0.3 to ~ 6 days, followed by a 
sharp drop at periods of > 8 days. Since the detection probabil¬ 
ity of BEER binaries decreases with orbital period (see Fig. [ 7 ]), 
this rise is probably real. In the past three decades, several stud¬ 
ies reported a log-normal period distribution of solar-like bina¬ 


ries, with a peak at lo gF ~ 5 (e.g.. lDuauennov & Mavoilll99ll : 
IRaghavan et akllToiOl) . and Fig. [TT] shows that our period his¬ 
togram fits well the log-normal fit of lRaghavan et al.l (1201 Ol) out 
to P ~ 8 days. The scaling factor between the two samples is 
1.61, which means that for a bin size of AlogP = 0.214, as we 
use here, we would obtain ~ 74 binaries at the peak of the dis¬ 
tribution at logP = 5.03. The histogram sharply drops at periods 
of > 8 days both because at such periods BEER approaches its 
sensitivity limit when applied to CoRoT long-run light curves 
and because the AAOmega observing run was limited to seven 
nights. 

The mass-ratio distribution of short-period binaries was the 
subje ct of several in-depth studies in the last three decades 


(e.gjHalbwac 


Duquennov & Mayor 


iHalbwachs et al.ll2003b [Goldberg e t al.ll2003t iFisher et al 


azeh et a! 


20051) . The main debates regarding the mass-ratio distribution of 
spectroscopic binaries are about the existence of peaks at q ~ 1 
and/or at q ~ 0.2 and about the shape of the distribution toward 
lower q values—whether it is monotonically increasing, decreas¬ 
ing, or flat. It is therefore of interest to plot the distribution of the 
assigned mass ratios of our sample, as was done in Fig. ITT1 
The histogram presents three features, which we describe 
from top to bottom. (1) A clear peak at q ~ 1. (2) The number 
of binaries increases with decreasing q. (3) The histogram peaks 
at q ~ 0.2 and then falls s harply towards q ~ 0.1. However, 
correcting our sample for the lQpikl (1 19241) effect might have sig¬ 
nificantly diminished th e q ~ 1 peak (e.g. JGoldberg et alll2003b 
IHalbwachs et al.ir2003l) . In addition, at q < 0.6 our h istogram 
is dominated by SBls, and iMazeh & Goldbergl (1 19921) have al¬ 
ready shown that using an average sin / value to all SBls, like 
we did here, might produce a monotonically increasing q distri¬ 
bution when the actual true distribution is flat. Lastly, the drop 
toward q ~ 0.1 can be explained both by BEER’s sensitivity 
for CoRoT long-run light curves, which probably approaches its 
limit at such low mass-ratios, and by the limited RV precision 
of our AAOmega data. Therefore, before reaching any definite 
conclusion, our data need an in-depth analysis to derive a real 
mass-ratio distribution of the CoRoT BEER sample (Shahaf et 
al., in prep.). 


8. Summary 

We have presented AAOmega RV confirmation of seventy new 
non-eclipsing short-period binaries found by BEER in CoRoT 
light curves. The confirmed binaries span a range of 0.3-10days 
in orbital period, showing a clear rise in the number of bi¬ 
naries per AlogP toward longer periods. Our spectral analysis 
shows that the primary stars in the detected binaries are typi¬ 
cally of spectral type G or earlier, and the mass ratio (includ¬ 
ing the mass ratio assigned to SBls) spans a range of 0.03-1. 
The mass-ratio histog ram resembles a double-peak distribution 
(e.g Halbwachs et al.l2003l) . but since we did not correct our re¬ 
sults, neither for some well-known selection effects (e.g., lOpikj 
1924 nor for the non-isotropic inclination distribution (e.g., 
Mazeh & Goldbergl 19921) . the underlying true distribution might 
also be flat. Nevertheless, the orbital solutions presented here can 
be used in the future for an in-depth study of the mass-ratio dis¬ 
tribution of short-period bin aries , similar to the study presented 
by iGoldberg. Mazeh. & Lathami ( 2003 /). 

On the lower end of our detection capability, we have de¬ 
tected two BD candidates on a ~ 1 day period orbit around 
F-G stars. We considered them as BD candidates and not bona 
fide BDs for two reasons: (1) the true inclination is not known, 
and (2) higher resolution and/or S/N spectroscopic observations 
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might reveal a faint companion that could not have been found 
in the AAOmega spectra. 

Relativi stic beaming was alrea dy detected in CoRoT light 
curves (e.g. jMazeh & Faiglerl20ldi) . but this is the first time that 
BEER detected non-eclipsing binaries in CoRoT light curves. 
To roughly estimate the expected number of such binaries that 
can potentially be detected with BEER in CoRoT long-run light 
curves we first need to correct for the 200 binaries that were 
selected for follow-up, but were not observed, taking into ac¬ 
count their BEER priorities, and then for the total number of 
CoRoT long-run light curves (~ 110,000), assuming that short- 
period binaries do not have a preference for center or anticenter 
fields. Doing so, we estimate that observing all 481 selected can¬ 
didates would have brought the number of confirmed binaries to 
~ 110, and that the expected number of beaming binaries that 
can be detected in CoRoT long-run light curv es is ~ 3 00. This 
confirms the prediction made by IZucker. Mazeh. & Alexander! 
(120071) : CoRoT and Kepler will be capable of detecting hundreds 
of non-eclipsing beaming binaries. 

Investigating BEER’s false-alarm probability and nature in 
CoRoT long-run light curves, we showed that red giants in¬ 
troduce a major source of false candidates and demonstrated a 
way to improve BEER’s performance in extracting higher fi¬ 
delity samples from future searches. Understanding the depen¬ 
dence of the fidelity and completeness of such future samples 
on parameters such as period, spectral type, and BEER-model 
score might enable using these samples to derive some statistical 
properties of the short-period binary population, like the period 
distribution, even before RV follow-up is done. Merging well- 
corrected large samples of CoRoT and Kepler beaming binaries 
with equivalent EB samples (e.g. JPrsa et al.|[20TTl) might create 
large and well-defined samples of short-period binaries, which 
may shed some light on binary formation and evolution. 

Radial-velocity follow-up observations of beaming and 
eclipsing binaries will continue to play a key role in these ef¬ 
forts because in most cases this is the only way to measure their 
mass ratio and the mass of the secondary companion. Scaling 
the AAOmega time that was required to make the detections 
reported here, we estimate that only three additional such runs 
would be required to detect virtually all beaming binaries that 
can be detected in CoRoT long-run light curves. This again 
demonstrates the efficiency of usin g mu lt i-obj e ct spec trographs 
for this task (e.g.. IZucker. Mazeh. & Alexanderll2007l) . One ad¬ 
vantage of the beaming binaries over the EBs is the ability of de¬ 
tecting binaries not only at edge-on inclinations, which widens 
the window for detecting intrinsically rare systems such as BD 
and massive-planetary companions to mai n-sequence stars (e.g., 
iHalbwachs et alJl2000biFaigler et al. 11201 3l) . 
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Fig. 12. Continued. 
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Fig. 12. Continued. 
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Fig. 13. AAOmega RVs and the best-fit Keplerian models (solid lines) of the confirmed BEER SB2s listed in Table[7] Primary RVs 
are marked with circles and secondary RVs with triangles. 
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Fig. 14. AAOmega RVs (black circles) and the best-fit Keplerian models (solid lines) of the variable components in the confirmed 
BEER diluted binaries listed in Table [8] 
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Fig. 15. Phase-folded and binned light curves (blue) and the best-fit BEER models assuming a circular orbit (green) of the 70 BEER 
binaries confirmed by AAOmega. 
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